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Abstract 



We study the two-photon process for the transitions «s — » Is and nd — > Is in hydrogen up to large n. For n < 20 we provide simple 
analytic fitting formulae to describe the non-resonant part of the two-photon emission profiles. Combining these with the analytic 
form of the cascade-term yields a simple and accurate description of the full two-photon decay spectrum, which only involves a sum 
over a few intermediate states. We demonstrate that the cascade term naturally leads to a nearly Lorentzian shape of the two-photon 
profiles in the vicinity of the resonances. However, due to quantum-electrodynamical corrections, the two-photon emission spectra 
deviate significantly from the Lorentzian shape in the very distant wings of the resonances. We investigate up to which distance the 
two-photon profiles are close to a Lorentzian and discuss the role of the interference term. We then analyze how the deviation of the 
two-photon profiles from the Lorentzian shape affects the dynamics of cosmological hydrogen recombination. Since in this context 
the escape of photons from the Lyman-a resonance plays a crucial role, we concentrate on the two-photon corrections in the vicinity 
of the Lyman-a line. Our computations show that the changes in the ionization history due to the additional two-photon process 

from high shell (n > 2) likely do not reach the percent-level. For conservative assumptions we find a correction hN^IN^ 0.4% at 

redshift z ~ 1 160. This is numerically similar to the result of another recent study; however, the physics leading to this conclusion is 
rather different. In particular, our calculations of the effective two-photon decay rates yield significantly different values, where the 
destructive interference of the resonant and non-resonant terms plays a crucial role in this context. We also show that the bulk of the 
corrections to the ionization history is only due to the 3s and 3d-states and that the higher states do not contribute significantly. 

Key words. Atomic processes ~ two-photon transitions - Cosmology ~ Cosmic Microwave Background 



1. Introduction 

During the epoch of cosmological hydrogen recombination (typ- 
ical redshifts 800 < z < 1600), any direct recombination of elec- 
trons to the ground state of hydrogen is immediately followed by 
the ionization of a neighboring neutral atom due to re-absorption 
of the newly released Lyman-continuum photon. In addition, be- 
cause of the enormous difference in the 2p Is dipole transi- 
tion rate and the Hubble expansion rate, photons emitted close to 
the center of the Lyman-a line scatter ~ 10*^ - 10^ times before 
they can finally escape further interaction with the medium and 
thereby permit a successful settling of electrons in the hydrogen 
ground state. It is due to these very peculiar circumstances that 
the 2s Is-two-photon decay process, which is ~ 10*^ orders of 
magnitude slower than the Lyman-a resonance transition, is able 
to substantially control the dynamics of cosmologica l hydrogen 
recombination dZeldovich et al.l [19681 |Peeblesl[l968h . allowing 
about 57% of all hy drogen atoms in the Unive rse to recombine 
at redshift z < 1400 (IChluba & Sunyaevl2006ah . 

The tremendous success in observations of the cosmic mi- 
crowave background t emperature and p olarization anisotropies 
(iHinshaw et all |2006'; 'Pag e et al.1 l2006h has recently motived 
several works on high precision computations of the cosmologi- 



2007at iHirata & Switzeij l2007t ISwitzer & Hiratal l2007bl: 
Kholupenko et al .I l2007h recombination history. 



cal hyd r ogen (iDubrovich & GrachevI 120051: IChluba & Sunvaeyl 
2006bt iKholupenko & IvanchikI l2006t iNovosvadlvil l2006t 
Rubino-Martfn et al.l 120061: IChluba et al l 120071: IWong & Scott 



2007t IChluba & Sunyaevll2007bl) and helium jSwitzer & Hirata 



One interesting additional phys ical process, which had b een 
neglected in earlier computations dSeager et al ] ll999ll200flh . IS 
connected to the two-photon transitions from high ns and nd- 
sta tes to the ground state of hy drogen and was first proposed 
by IDubrovich & GrachevI (l2005h . In their computations a simple 
scaling for the total two-photon decay rate of the s and d-states 
in hydrogen was given and, including these additional channels 
leading to the Is-level, corrections to the ionization history were 
found that exceed the percent-level. These modifications would 
have a strong impact on the determination of the key cosmolog- 
ical parameters (Lewis et al. 2006) and therefore require careful 
consideration. 

More recently, theoretical values for the non-resonan^ two- 
ph oton decay ra t es of the 3s and 3d level based on the work 
of ICresse r et al. (1986^ were u tilized to improve the computa- 
tions of [Dubrovich & GrachevI (12005.) . showing that the effect 
of two-photon transitions from higher levels o n the recombi- 
nation history is likely to be less than ~ 0.4% dWong & ScottI 
120071) . However, to our knowledge, neither extensive calcula- 
tions of the two-photon decay rates for the transitions ns — » 
Is and nd — > Is from high levels exist nor detailed re- 
ports o f direct measur e ments can be found in the literature, 
so that IWong & ScottI ( |2007|) also had to extrapolate exist- 
ing values towards higher levels, largely relyin g on the pre- 
vious estimates by IDubrovich & GrachevI d2005l) and the val- 
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' We base our definition of this terminus on the energies of the in- 
volved intermediate states (see Sect.[2]for details). 
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ues given in ICresser et al. I (119861) . Here we argue that, due 
to quantum-interference, it is difficult to separate the con- 
tributions of the pure two-photon process from the reso- 
nant single photon processes. Therefore to answer the ques- 
tion how much two-photon processes are affecting the re- 
combination history requires a more rigorous treatment in 
connection with the radiative tr a nsfer and escape of photons 
('Varshalovich & Svunvaev 1968; Grachev & Dubrovich '1991^, 
[R ybicki & dell' Antonio 1994 ; Chluba & Sunyaev 2007 a) from 
the main resonances (especially the Lyman-a Une). Also 
iHirata & Switzej ( |2007|) re-analyzed the importance of the two- 
photon process in the context of cosmological helium recombi- 
na tion and showed that, f or high values of n, the rate estimates 
bv lDubrovich & Grachev! (120051) are rather rough and that in par- 
ticular the linear scaling with n fails. Here we provide some con- 
servative lower limits on the possible impact of the two-photon 
transitions on the hydrogen recombination history. We show that 
the non-resonant contribution to the two-photon decay rate in- 
deed scales oc n (see Sect. 14.3.11 ). However, due to destructive 
interference between the resonant and non-resonant terms, the 
effective two-photon decay rate is much lower and actually de- 
creases with n. 

If one considers an isolated neutral hydrogen atorrQ with the 
electron in some excited state («, /), then because of the finite 
lifetime of the level, the electron will reach the ground-state after 
some short time (typically ~ 10"'^ s), in general releasing more 
than one photon. Astrophysicists usually describe this multi- 
photon cascade, also known as Seaton-cascade, as a sequ ence of 
independent, single-step, one-photon processes (e.g. see ISeatonI 
Il959h . where every resonance has a pure Lorentzian shape. This 
approximation should be especially good in the presence of 
many perturbing particles (free electrons and ions), such as in 
stellar atmospheres, which destroy the coherence of processed 
involving more than one intermediate transition. However, in ex- 
tremely low density environments, like the expanding Universe 
during cosmological hydrogen recombination , hardly any per- 
turbing particle is with in the Weisskopf -radius (IWeisskopfll932t 
ISobelman et al.in"995b . so that the coherence of two-photon and 
possibly multi-photon transitions is maintained at least for the 
lower shells. Here we consider the simplest extension to the clas- 
sical treatment of the multi-photon cascade and focus only on the 
two-photon process. 

Beginning with the paper of iGoppert-Maverl ( 193 1|). several 
textbooks of quantum electrod ynamics (lAkhiezer & Berestetskiil 
119651; iBerestetskii et al.|[T982l) discuss the two-photon emission 
process. Using the formulation of quantum-electrodynamics, 
one naturally obtains a nearly Lorentzian shape of the line pro- 
files in the vicinity of the resonances, which also allows us to 
check for tiny deviations in the real two-photon emission from 
the spectrum obtained using the simplest classical cascade treat- 
ment. As we discuss below, at least in the decay of high s and d- 
states quantum-electrodynamical corrections lead to additional 
broad continuum emission and strong deviations of the profiles 
from the natural Lorentzian shape in the very distant wings of the 
resonant lines. In this paper we investigate up to what distance 
the wings of the two-photon emission spectrum in the vicinity 
of the Lyman-ff line continue to have a Lorentzian shape. These 
deviations in the red wings are the reason for the corrections to 



- We restrict ourselves to the non-relativistic formulation of the hy- 
drogen atom and assume that it is at rest in the lab frame. 

^ More rigorously, collisions more likely destroy the coherence for 
photons emitted very close to the line center, since in the distant wings 
the typical characteristic time of pro cesses involving the correspon ding 
quantum-states is much shorter (e.g. lKarshenboim & Ivanovll2007h . 



the hydrogen recombination history due to the two-photon tran- 
sitions from high s and d-states. Similarly, these modifications 
of the Lyman-ff line profile could also be important during the 
initial stages of reionization in the low-z Universe. One should 
mention that the developed picture is valid for the primordial 
chemical composition of the Universe, which is characterized by 
a complete absence of heavy elements, e.g. dust and low ioniza- 
tion energy species that would influence the escape of Lyman- 
a photons in planetary nebulae and H ii regions in present-day 
galaxies. 

For hydrogen, several publications on the theoretical value 
of the total 2s —> Is two-photon decay rate can be found 
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1986; Goldmanlll989l) with recent computations performed bv 



Labzow skv et al J (!2005l) yielding Aasis = 8.2206 s . In these 
calculations one has to consider all the possible intermediate 
states (bound and continuum) with angular momentum quan- 
tum number 1-1, i.e. p-states. Within the non-relativistic treat- 
ment of the hydrogen atom for the metastable 2s-level, no p- 
state with energy lower than the 2s-state exists; hence the two- 
photon process only involves transitions via virtual intermediate 
states, without any resonant contributions. Therefore the total 
two-photon decay rate of the 2s-level is very low and the 2s- 
state has an extremely lo ng lifetime (~ 0.12 sec). We show that 
the formulae obtained bv !Cresser et al.! ( !l986!) are not applicable 
in this case (see Sect. 14. II ). 

Also, some calculations for the 3s and 3d t wo-photon transi- 
tions t o the ground state h ave been carried out (!Ouattropani et al.! 
!1982!; !Tung et al.! ! 19841; !Florescul !l984l) . but here a problem 
arises in connection with the contribution from the intermediate 
2p-state, which has an energy below the initial level. The corre- 
sponding term is dominating the total two-photon decay prob- 
ability for the 3s and 3d two-photon process and is connected 
with the resonant transition via an energetically lower level. It 
can be interpreted as a cascade involving the quasi-simultaneous 
emission of two-photons. As in the case of the 2s-level, a broad 
continuum emission also appears due to transitions via virtual 
intermediate states, with energies above the initial level. This 
continuum is not connected to any resonances and therefore has 
a much lower amplitude. In addition to the cascade-term and this 
non-resonant term, an inte rfe rence-teim also appears for which 
a clear interpretation is difficult within the classical formulation. 
Similarly, in the two-photon decay process of higher ns and nd- 
states to the ground state, (2« - 4)-resonances appear, yielding 
complex structures in the distribution of emitted photons. Some 
addition al examples of 2 y er nission spec t ra caii also be found in 
Ouattropani et alj ( !l982l) and !Tunget all ( !l984 . 

As mentioned above, astrophysicists usually interpret the 
two-photon cascade as a 1 + 1 -single photon process. Since even 
in the full two-photon formulation, the cascade-term dominates 
the total two-photon decay rate (hence defining the lifetime of 
the initial «s and «d-states), in a vacuum the total two-photon 
decay rate should be very close to the 1 -i- 1 -single photon rate 
of the considered level. In the 1 + 1 -photon picture, the sponta- 
neous two-photon decay rate is simply given by the sum of all 
spontaneous one-photon decay rates from the initial state, since 
after the detection of one photon, say a Balmer-a photon in the 
3s— >ls transition, in vacuum the presence of a Lyman-a pho- 
ton is certain and therefore should n ot affect t he tota l 3s decay 
probability. For the 3s and 3d-states !Florescul (! 19841) computed 
the total spontaneous two-photon decay rate and indeed found 



6.469 X 10^ S-' 



6.317 X lOS- 



and A~l , 

3d— >ls 



,p . . . „ j(d->ls 3d->2p 

This also suggests that, very close to the reso- 
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nances, the 1 + 1 -photon description provides a viable approx- 
imation, in which the line profile is very close to a Lorentzian. 
However, as we show below, quantum-electrodynamical correc- 
tions (e.g. virtual intermediate states, interference, correlations 
of the photons in energy) lead to differences in the two-photon 
profiles in comparison with the 1 + 1 -single photon profile, which 
are significant especially in the distant wings, far from the reso- 
nances. In particular the interference term plays a crucial role in 
this context and cannot be neglected. 

2. Two-photon transitions of the m and nd-states 

Considering only cases when the initial states is either an «s 
or «d-level and the final state corresponds to a s-level, one can 
simplify the general formula for the t wo-ph oton transition prob- 
ability as given by iGoppert-Maverl (Il93lb considerably. First, 
the average over the directions and polarizations of the emitted 
photons can be carried out immediately, since within the non- 
relativistic formulation, one can separate the radial and angu- 
lar parts of the wave function. For s — > s-transition, this leads 
to a global factor of 1 /27, w hile this average yields 2/135 for 
d — > s-transitions (see lTung e t al. 1984i). Afterwards, the proba- 
bility for the decay «,/,■ — > Is (where /,■ = or /, = 2) with the 
emission of two photons can be written in terms of the integrals 
{R„'i'\r\Rni) = {R„'i'\r\R„i)* = 7?,,-/- r^7?„/ dr over the normal- 
ized radial functions, R„i{r), for whic h expUcit expressions ca n 
be found in the literature (e.g. §52 in iBerestetskii et aTlll982h . 
Then the probability of emitting one photon at frequency v and 
another at v' in the transition (n,7,) — > Is is given by 



(la) 



M^J] {Ru\ r \R„p) {r„p\ r \R„,i) g„{v) (lb) 



gn(y) 



1 1 

+ 



El, - E„ + hv E„ - E„ + hv' 



(Ic) 



Here C; = 9q'*'c7?h/2'" x fl//v^^;^ ~ 4.3663 s-' x fl//v^^j^, with 
ao = 1 and ao—llS, and where a is the fine structure constant, c 
the speed of light, V2sis the 2s-ls transition frequency, and R^ the 
Rydberg constant for hydrogen. Due to energy conservation, the 
frequencies of the emitted photons are related by v -n v' = v/is, 
where here v,is is the transition frequency between the initial 
level and final Is-state. The energy of level « is given by 
E„ - -Eis/n^, with the Is-ionization energy of the hydrogen 
atom Els a: 13.6 eV. The sum in Eq. ([T]) has to be extended by an 
integral over the continuum states, but for convenience we omit 
these in our notation. 

For «, > 2 and n < «,, it is clear from Eq. ([T]i that at 
hv = £■„, - En and hv = E„ - Ei^, i.e. corresponding to the 
resonance frequencies to energetically lower levels, one of the 
denominators inside the sum vanishes, leading to a divergence 
of the expression. As we discuss below (Sect. 12.1.3b . including 
the lifetime of the intermedia te states provides a possibility of 
removing these singularities (lLowlll952l: iLabzowskv & ShoninI 
12004 : however, a consistent consideration of this problem re- 
quires a more sophisticated treatment beyond the scope of this 
paper Physically, transitions to intermediate states with energies 
En > are virtual We split up the sum over all the inter- 
mediate states where - IjE„>e„ and 



Strictly speaking, photons appearing in the distant wings 
of the resonance are also connected with virtual states (e.g. 
iKarshenboim & Ivanovll2007h . 

' Within the non-relativistic treatment of the hydrogen atom, this is 
equivalent to separating the cases n < n, and n > Since the energy 



2i = YjE„<e„ denote the sum over virtual and real intermediate 
states, respectively. Then we can write 



(2) 



Henceforth, we refer to the contribution to the transition ma- 
trix element from as the cascade part, due to l2vP as non- 
resonant pait, and the rest as the interferenc^parX. This distinc- 
tion is ad hoc and not unique, but only motivated by our separa- 
tion of the infinite sum, which will turn out to be fairly conve- 
nient in terms of evaluation of the spectrum and total two-photon 
decay rate. 

The non-resonant contribution to the two-photon decay 
probability, dW™, is then given by 

dC;/,^is=C,, ^'^"l^mf dv (3a) 

OO 



For the cascade and interference terms, one instead has 



(4a) 



Meas = ^ (/?lsk \Rnp) (^«pl r \R„I,) g„iv) (4b) 



and 



n=2 



^Kl^l^ = y'y" Mnr + M,,, Ml] dv, (5) 



respectively. 



2.1. Total two-photon decay rate 

In order to obtain the total two-photon decay rate in vacuum one 
now has to integrate Eq. ([T]i over all possible frequencies v. The 
corresponding integral can be cast into the form 



iy)dy, (6) 



with y = v/v,is and where the factor of 1/2 avoids double- 
counting of photons. In Eq. Q we have introduced the two- 
photon decay profile function, 

<f>%,^H(y) - €.,^u(y) + C-is(3') + CUis(3')' (7) 

which is the sum of the two-photon profiles due to the non- 
resonant, cascade, and interference terms, respectively. Below 
we discuss the contribution of each term separately. 



2.1 .1 . Interpretation of the two-photon emission profile 



Physically the two-photon emission profile or spectrum (p^^j.^i^ 
defines the number of photons that are released per second in 
the frequency interval between v and v + dv. If one integrates 
over the whole spectrum, this therefore yields the total number 



of the continuum states is always greater than for the bound states, the 
former are associated with the case n > n,. 

* In principle one should be more accurate by calling this contribu- 
tion resonant/non-resonant-interference term, since also some level of 
interference is already included inside the resonant and non-resonant- 
term, which is absent in the 1 + 1 -photon picture. However, we generally 
do not make this distinction. 
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of photons emitted per second due to the two-photon transition. 
The two-photon profile includes both photons at the same time, 
so that the total two-photon transition rate per initial s or d-state 
has to be divided by a factor of 2 (see Eq.|6]l. 

Due to energy conservation it is clear that, when detecting 
a photon that was produced in a particular two-photon transi- 
tion from some initial «s or «d-state to the ground state, at a fre- 
quency V, the other photon has frequency v' = v„is - v. Therefore 
also the probability to release a photon at v should be equal to 
the probability for the emission of a photon at v„is - v, a prop- 
erty that is reflected in the symmetry of the two-photon profiles 
around y =1/2 (see Sect. [3]for more explicit examples). 

2.1 .2. Total two-photon decay rate due to the non-resonant 
contribution 

In Eq. O we have introduced the non-resonant two-photon de- 
cay profile function, ip™i^iJj), which can be written as 



CO 

M„ = 2 (^Isk \Rnp) {R„p\ r \R„,,) fniy) 

1 1 

fniy)-- 



y+yT y-yT 



(8a) 
(8b) 

(8c) 



with G,ij/; — vfjj,C/;. Here we defined the frequencies 



yr - 



n\nj- 1) 



nHn]- 1)' 



(8d) 
(8e) 



which for n > ii/ vary within the ranges < y™'^ < l/{nj - I) 
and 1 < >';f'+ < «2/(«2 - 1). For n = n,- one finds = 
and y™'^ - 1 . Since the sum in Mm only involves intermediate 
states with « > «,, the denominators of /„ never vanish within 
the interval < y < 1, and for « - iii the factors y^ and (1 - y)^ 
ensure that Mm approaches zero within the limits y — > and 
y — » 1 . In addition is real and symmetric around y = 1/2. 

To compute the total rate one now only has to replace 
(p^^^i^iy) in Eq. Q, by the expression To evaluate the sum 
and integrals over the radial functions, we used Mathematica. 
Normally we restrict ourselves to the first 200 terms in the sum, 
but computations with up to 4000 terms were also performed for 
the 2s, 3s, and 3d rates. Within the assumptions the results for 
the other levels should be correct to better than ~ 1%. To make 
cross checks easier, we give the expression for the necessary ra- 
dial integrals {R„'i' \ r \Rni) up to «, = 5 in AppendixlAl 

2.1 .3. Total two-photon decay rate due to the cascade and 
interference terms 

For the cascade and interference terms special care has to be 
taken close to and at intermediate distances from the resonance 
frequencies hv - E„. - E„ and hv = E„ - £15. As mentioned 
above, a consistent treatment of this problem requires more so- 
phisticated methods, including the amplitudes of sev eral addi- 
tional processes (e.g. Karshenb oim & Ivanov, l2007l) . than are 
within the scope of this paper One simple approximate solu- 
tion to this problem can be given when taking the lifetime of 



the intermediate states into account as a small imaginary con- 
tribution to their energjQ. Including this shift into the equa- 
tions for the 2p-transition leads to the classical expression of the 
Lorentzian within the formulation of Quantum Electrodynamics 
and can be attributed to the first order radiative corrections of th e 
one-photon pro cess dLowl Il952t iLabzowsky & ShoninI l2004l) . 
iFlorescul (1 19841) used this approach to compute the total 3s and 
3d two-photon decay rates and simply replaced the energy, E2p, 
of the 2p-state by, £'2p - ih^ip/'^, where is the width of 
the 2p-state due to spontaneous transitions. Except for the 2s- 
state, summing all the one-photon dec ay rates (e.g. these values 
can be computed using the routines of lStorey & HumrneJll99lh 
should yield a very good approximation for the total lifetime of 
any given initial level in the hydrogen atom. For estimates, we 
therefore follow this approximate procedure and replace the en- 
ergies of all intermediate p-states by Ei,p — > E^j, - iliTlp/2. Here 



^ np — ^ 

ate «p-state due to spontaneous transitions. Note that for n > 2 
the p-states can decay via channels, which do not directly lead 
to the Is-state, thereby leading to the possible emission of more 
than two photons. 

The two-photon contribution to the total lifetime of the s and 
d-levels for transitions to the ground state should be close to the 
value following from the sum of the rates for all possible ly- 
transitions to lower-lying intermediate p-states multiplied by the 
probability that the electron will 'afterwards' go directly to the 
Is-leveS 



„„ _ ,i,,„ = Y,,,,, a''' ,„ is the total ly-width of the intermedi- 



«s/d 



n-1 

n'=2 



«s/d^«'p 



n'p 



(9) 



where ^j,,/;,^,,/,, is the spontaneous one-photon Einstein coef- 
ficients of the considered transition from the initial s/d-state0. 
Although for the s-levels, this already accounts for all the possi- 
ble one-photon decay channels, d-states with n > 4 can also de- 
cay via intermediate f-states. However, in this case again more 
than two-photons have to be released in order to reach the Is- 
level and therefore do not contribute anything here. The addi- 
tion (|9]l to the energies is not important for the non-resonant 
term, but it becomes crucial very close to the resonances for 
both the cascade and interference term. However, here we are 
mainly interested in the behavior at large distances from the 
poles (|Av| » r',p), so that this replacement has no direct effect 
on the discussed cosmological application (see Sect. 15.3b . 

Here one may ask, why the lifetime of the initial level is not 
includedl Physically thi s is motivated by the idea tha t, follow- 
ing the interpretation of IWeisskopf & WigneJ (Il930l) . we con- 
sider one particular initial 'energy sub-level' and do not specify 
the process that populated it. Therefore the final profile should 
be independent of the shape of the distribution of energy-sub- 
levels around the mean energy of the initial state. One can also 
consider this as equivalent to neglecting any possible reshuffling 
of the electron by perturbing particles while it is in the initial 
state. However, in the computation presented below we do not 
approach the resonances so close that these difl'erences would 
play any role. 



' We neglect the small correction to the real part of the energy caused 
by the Lamb-shift. 

The authors wish to thank E.E. Kholupenko for pointing out some 
inaccuracies related to this aspect. However, the modification did not 
affect the results of this paper. 



Here A 



ly 



refers to either the A''' ,„ or the A'*! ,„ decay rate. 



ly 



ns/d— >«'/' 

We use this notation more often below, 
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Figure 1. Non-resonant two-photon emission spectra, Eq. ([8]), 
for several transitions. All curves are normalized to unity aty - 
1/2. The values of 0"(l/2) can be found in Appendix IbI 



With the notation of Sect. 12.1.21 we now introduce the cas- 
cade and interference two-photon emission profiles by 



where the cascade matrix element is given by 

«/-! 

Meas = {Rh\ r \R„p) {R„j>\r\R„,,) K(y) 



n=2 



Kiy) = 



1 



1 



(lOa) 
(10b) 

(10c) 
(lOd) 



y - y'n"" ■ - iS„ y-yn ' - iS„ 

Here 6„ = r„p/4;rv,is accounts for the energy-shifts due to the 
finite lifetime of the intermediate «p-state. We have also intro- 
duced the resonance frequencies 

(11a) 
(lib) 

Since for the cascade and interference term n < n,, these now 
have values strictly within the range Q < y < I. 
Defining the function 

a 



yT'-^y^-- 



L{a, b) = 



(12) 



which actually has the generic shape of a Lorentzian, then 
fniy) = r„(y) + ifniy), where the real and imaginary part of /„ 
are defined by 

/^(y) = L(y - f:'-\ 6n) - L(y - y^'^ S„) (13a) 
fl(y) = L{6„,y- y^'^) - L{S,„y- yT''), (13b) 
respectively. Introducing k„ = (^Ri^\ r \R„^ {Rnp\ r l^^n,/,) one can 
rewrite |McasP and |Mi„tP = M^AK^s + ^cas] as 

n,-l «-l 



;i=2 



n=2 m=2 



"/-I 



«=2 



(14b) 



To compute the total rate one now only has to replace 'p'^Ji^^J.y) 
in Eq. (|6]l, by the corresponding expressions ( fTOl i. 



3. Two-photon emission spectra 

3. 1 . Non-resonant two-photon emission spectra 

In Fig. [U we present the profile functions for the non-resonant 
contribution to the two-photon decay spectrum. All the pro- 
files have a maximum at y = 1/2. For the «s — > Is-emission 
profiles the difference in the shape of the curves is quite big, 
while for initial d-states in general the profile does not vary as 
much. However, in both cases the amplitude at y = 1/2 changes 
strongly, increasing towards larger n (see AppendixlBl). 

Due to our separation of the infinite sum over the intermedi- 
ate substates, the sums in the cascade and interference terms be- 
come ^n;fe. This allows us to evaluate Mnr numerically and use 
convenient fitting formulae for their representation. Realizing 
that is symmetric around y = 1/2 and that it scales like 
~ Ijy and ~ 1/(1 - y) at the boundaries, we approximated 
y{ 1 -y) Mni . In AppendixlBl the obtained formulae for all «s ^ Is 
and «d — > Is transitions up to « — 20 are given. For the non- 
resonant term within the range 10"^^ < y < 0.999, these approx- 
imations should be accurate to better than 0.1%. Since these are 
fast and simple to evaluate they should be useful for analytic es- 
timates and numerical applications. Note that all non-resonant 
matrix-elements are negative. 



3.2. Tine 3s and 3d two-pinoton emission spectrum 

In Fig. 12] we give the full two-photon emission spectra for the 
3s — » Is and 3d — > Is transition. In both cases, the non- 
resonance term increases the wings of the profiles close to y ~ 
and y ~ I. However, the interference between the non-resonant 
and cascade part is destructive in the central region {y ~ 1/2) and 
significantly reduces the amplitude of the total two-photon emis- 
sion. It even leads to full cancellation at y ~ 0.22 and y ~ 0.78 
for the 3s-level, whereas for the 3d-level the photon production 
does not vanish in the region between the resonances (see also 
iTung et al.lfT984l) . 

For both the 3s —> Is and 3d ^ Is transitions, only one 
term in the cascade appears, which is related to the transition 
via the intermediate 2p-state. The matrix elements for these are 
A-^^ ^ 1 .21 1 and k^'^ » 6.126, and according to Eq. (fTTT i the reso- 
nance frequencies are at yj"^'^ = 5/32 (Balmer-o- transition) and 
y^"^'' - 27/32 (Lyman-a transition). With the Eqn. given in Sect. 
12.1.31 and using the fitting formulae according to Appendix |B] 
one can analytically approximate the full two-photon emission 
spectrum. As Fig.|2]shows the agreement is excellent at all con- 
sidered frequencies. 

3.3. The ns and nd two-photon emission spectrum 

As an example, in Fig. [3] we present the fuU two-photon emis- 
sion spectra for the 5s —> Is and 5d ^ Is transition. Again 
one can see that the interference term strongly affects the shape 
of the spectrum in the wings of the resonances. In particular, 
destructive interference close to y = 1/2 strongly reduces the 
total amplitude of the emission. For the 5s-level, interference 
leads to full cancellation of the photon production (y ~ 0.28 
and y ~ 0.72) in the region between the innermost resonances, 
whereas the photon production does not vanish within this range 
for the 5d-level. This difference is characteristic of the shape of 
the s and d-two-photon spectra, also for higher values of n. It is 
also clear, that using the Eqn. given in Sect. 12.1.31 together with 
the fitting formulae according to Appendix |B] one can analyti- 
cally approximate the full two-photon emission spectrum with 
very high accuracy in the full range of considered frequencies. 
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Figure 2. Two-photon emission spectra for the 3s ^ Is and 3d ^ Is transitions. The non-resonant, cascade, and combined spectra 
are shown as labeled. Also we give the analytic approximation based on the equations discussed in Sect. 12.1.31 together with the 
fitting formula for the 2s and 3d non-resonant matrix element (see Appendix|B]). The resonances correspond to the Balmer-or and 
Lyman-a transitions. 



For initial states with a higher value of n,, more resonances 
(in total 2 n, -4) appear, but otherwise the spectra look very sim- 
ilar and do not add any deeper physical aspects. We checked the 
analytic approximations for the full two-photon emission spec- 
trum of several ns — » Isandnd — » Is two-photon transition up 
to n - 20 and always found excellent agreement with the results 
from our full numerical treatment. 



4. Total two-photon decay rates 

4.1. The two-photon rate for the Is-state 

It is clear that, within the non-relativistic formulation for the 
2s^ls-two-photon transition of hydrogen-like ions, only non- 
resonant contributions to the total lifetime exist. When we use 
Eq. Q together with Eq. ([S) and include the first 4000 terms 
in the infinite sum, we obtain the value A2sis = 8.2293s"', 
which fully agrees with the res ult of earlier comput ations 
(iBreit & Telle Al 940'; 'Kippei^'1950'; 'Spitzer & Greensteir]l[ " 

Klarsfeld 1969; John son 1972; Goldm an & D rake 1981;! 

1986HGoldmanll989l:lLabzowskv et alJ2065h . With the approx- 
imate formula for the non-resonant two-photon emission spec- 
trum for the 2s-state as given in Appendix|B] we obtain A2sis - 
8.2297s"', which shows the high accuracy of the approximation. 

Equation (|6]l, togeth er with Eg. (I p, is very similar to 
Eq. (14) in the work o f ICresser et al.l ( Il986h . The expression 



^ q- 

of ICresser et al 1 (Il986h was obtained using general arguments 
about the total lifetime of the considered level, and according to 
their work it should be applicable to all s and d-states of hydro- 
gen, yielding the two-photon correction to the lifetime. However, 
if applied to the hydrogen 2s-level one finds Ajjj^ = 1.4607s"' 
instead of A2sis - 8.2293s"'. The difference stems from the fact 
that here, like in the publications mentioned above, we included 
the term with n = tii in the sum Eq. (O. This shows that the 
largest contribution to the total 2s-two-photon decay rate (in this 
case equivalent to the non-resonant contribution) actually comes 
from the transition via the intermediate 2p-state, i.e. the matrix 
element ^.Risl I^2p) (^2pl I^2s) ~ -6.704, and cannot be ne- 



glected. This suggests that the arguments by lCresser et al ] (ll986h 
are incomplete, or at least not generally applicable. 

Although to our knowledge only rough direct measure- 
ments of the two-photon decay rate exist for the hydrogen 
2s-state (Kruger&Oed 1975; Cesa r etal . 199^, one can find 
experimental con firmations (IPriod [19721: iKocher et al.l 119721: 
iHinds et al.l [T978I) of the theoretical value for the two-photon 
decay rate of the hydrogen-like helium ion (Fth ~ 8.229Z^s"'), 
which do reach percent-level accu racy. Also m easurem ent for 
hydrogen-like Ar, F, and O exist ( Marrus & Schmieden 119721: 
ICocke et al.lll974HGould & MarruslTTgslT but with lower accu- 
racy. These experimental confirmations further support the idea 
that, in theoretical computations of the total two-photon decay 
rate and in particular the correction to the one-photon lifetime, it 
is not enough to consider only intermediate states with energies 
E„ > Ei. 

For hydrogen-like ions, care should be taken when comput- 
ing the 2s-two photon decay rate within the relativistic treat- 
ment. In this case the 2pi/2 level due to the Lamb-shift and 
fine-structure splitting energetically lies below the 2s 1/2 level. 
Increasing Z will make this shift even bigger, but as the mea- 
surements for He and Ar show, this intermediate state cannot 
contribute beyond the percent-level to the total lifetime of the 
corresponding 2si/2-state. This is also expected because the life- 
time of the 2s-state should not be strongly altered by the slow 
2si/2 — » 2p 1/2 transition (~ 1.6 x 10"^ s"'). In addition, the poles 
due to this intermediate state lie very close to v — > and v v,is 
and are therefore suppressed by the factors of v^'v'^ in Eq. ([T]i- 

4.2. The two-photon rates for the 3s and 3d-states 
4.2.1. The non-resonant contribution 



Using the formula given by ICresser et alJ (Il986l) . i.e. explic- 
itly neglecting the transition via the intermediate 3p-state, 
we can reproduce their values for the non-resonant contribu- 
tion to the two-photo n decay ra t es of t he 3s^ls and 3d— >ls 
transitions. Later Florescu et alj ( 1 19881) computed the se val- 
ues again within the framework of Cresser et al.l ( Il986l) but to 
higher accuracy. We are also able to reproduce these results 
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Figures. Two-photon emission spectra for the 5s ^ Is and 5d — > Is transitions. The non-resonant, cascade, and combined spectra 
are shown as labeled. Also we give the analytic approximation based on the equations discussed in Sect. 12.1.31 together with the 
fitting formula for the 5s and 5d non-resonant matrix element (see Appendix IbTi. The resonances correspond to the Brackett-a, 
Paschen-jS, and Balmer-y transitions at y < 1 /2, and Lyman-o', Lyman-/?, and Lyman-y aty > 1/2. 



= 8.22581 s-i and A^Ij'^'j^ = 0.131814 s"') up to a// 
given figures. Although the discussion in the previous sections 
has already shown that these values probably have no direct re- 
lation to the total corrections in the lifetime of the level due to 
the two-photon process, we computed them to check our own 
computational procedure. 

However, returning to our definition of the non-resonant two- 
photon decay rate, the transition via the intermediate 3p-state 
has to be included. We then obtain A™^^^^ = 10.556(10.558) s"' 
and Ajjj^j^ = 7.1474(7.1475) s"', where the values in paren- 
thesis were computed by integrating our analytic approxima- 
tion. In particular, for the 3d-level, this increases the non- 
resonant contribution to the total two-photon decay rate by a 
factor of ~ 54. If in the sum dSbt we only consider the term 

(^isk|/?3p)(^3pkl^3d)/3 * -5.199/3, with the function /3 = 
y"' -i-(l -y)"' and the integral y^(l -y)^/^ dy - 1 /6, then one 
obtains A?". , « 9.199s This shows that indeed the main 

3d— >ls 

contribution to the non-resonant part of 3d-two-photon decay 
rate comes from the transition via the intermediate 3p-state. 



4.2.2. The cascade and interference terms 

With the formulae given in Sect.|2l it should be possible to com- 
pute the total two-photon decay rate of the 3s and 3d-states. 
Including the lifetime of the intermediate 2p-state as discussed 
in Sect. 12.1.31 we also computed th e total lifetime of the 3s 
and 3d-states, and, in agreement with lFloresc ul (11984 . obtained 
values that were very close to the one expected from the one- 
photon lifetime. But as mentioned in Sect. 12.1.31 within the sim- 
ple approximation used to regularize the cascade and interfer- 
ence terms, it is not possible to compute the total correction to 
the one-photon lifetime, consistent in the considered order of the 
fine-structure constant a. In addition, as we will see in Sect. |5] 
this is not necessary for our cosmological application. 

However, in order to compare with other computations, it 
may be useful to give some additional intermediate results. We 
therefore also integrated the contribution of the interference 
term separately, yielding A'3"'^ J ^ = -10.810(-10.810)s"' and 
" -30.019(-30.019)s-i. This shows that, because of in- 
terference, the small increase of the decay-rate due to the non- 
resonant term (see Table [U is completely canceled, again em- 
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Table 1. The non-resonant contribution to the total two-photon 
rates for the transitions n,s — > Is and «,d — » Is up to «, = 20. 



We find that for «, > 20 one can use 



n, 








n,s->ls 


A"', , 

»ls 


2 


8.2293 s-' 


- 


- 




- 


3 


10.556 s-i 


7.1474 s-i 


12 


34.096 s-' 


29.687 S-' 


4 


11.951 S-' 


11.942 S-' 


13 


37.552 S-' 


31.410 s-' 


5 


13.741 S-' 


15.331 S-' 


14 


41.076 S-' 


33.114 s-' 


6 


15.954 s"' 


18.004 s"' 


15 


44.659 S-' 


34.801 S-' 


7 


18.501 S-' 


20.293 S-' 


16 


48.290 S-' 


36.478 s-' 


8 


21.301 S-' 


22.362 S-' 


17 


51.964 s-' 


38.146s-' 


9 


24.296 s-i 


24.296 S-' 


18 


55.674 S-' 


39.806 s-' 


10 


27.448 S-' 


26.144 s-' 


19 


59.416 s-' 


41.462 S-' 


11 


30.722 S-' 


27.935 s-' 


20 


63.185 S-' 


43.113s-' 




Figure 4. Non-resonant contribution to the total two-photon de- 
cay rate in vacuum for the ns and «d-states of the hydrogen atom. 
The results were computed using the first 200 terms above 



phasizing how important the interference term is. In Table[T]we 
included a maximal number of summands above the initial state, 
which was «sum = 4000 for «, € {2, 3) and «sum = 200 for «,■ > 3. 



4.3. The two-photon rates for the ns and nd-states 

For future computations and more complete considerations of 
the higher order correction to the lifetime of the «s and «d-states, 
here we now give the results for the total contribution of the non- 
resonant term to the two-photon decay rate. This contribution 
does not depend on the treatment of the poles in the cascade and 
interference terms. However, these values should have no direct 
relation to the total two-photon correction of the lifetime, but are 
mainly meant for cross-checks. 



4.3.1. The non-resonant contribution 

In Table [1] we summarize the values for the non-resonant con- 
tribution to the two-photon rates for the «s and «d-states up to 
«, = 20. The dependence of the non-resonant contribution to 
the total two-photon decay rate on «, is presented in Fig.H] For 
large n,- in both cases, the rates scale roughly linear, increasing 
towards larger n,. The slope is slightly steeper for the s-states. 



Am 



-15.857 s"' + 3.930 s"'«,- (15a) 
10.432 s"' + 1.636 s"'«/ (15b) 



within percent accuracy up to «, ~ 40. Explicitly computing the 
values for «, = 50 («sum = 200), we find A^;,^^,^ = 181.74 s-' 
and A"qj^j^ - 92.19 s-', using our full numerical treatment, 
and 180.64 s-' and 92.23 s-' with the approximations (fTsT l. We 
did not check up to which value of «, the formulae ( fTSl ) are 
applicable. Also, one should bear in mind that, above some 
value of n, » 1, the usual dipole approximation for the transi- 
tion matrix elements br eaks down jPubrovich & Grachevll2005l : 
iHirata & Switzeill2007h and other methods should be used. 

The linear scaling of the non-resonant contrib ution to the 
two-photon decay rate for «, » 1 was expected dPubrovichl 
Il987t iDubrovich & Grachevll2005h . but here we have included 
all virtual intermediate states in the sum. However, one should 
keep in mind that, due to the interference term, it is difficult to 
interpret this contribution separately. 



5. Astrophysical application 

5. 1 . Two-photon process in the single photon picture 

As described in the introduction, the standard procedure for 
treating the atomic transitions of electrons involving more than 
one photon is to break them down into independent, single-step, 
one-photon processes. This approximation should be especially 
good in the presence of many perturbing particles (free electrons 
and ions), such as in stellar atmospheres, which destroy the co- 
herence of processes involving more than one transition. Here 
we now explain how the two-photon process can be formulated 
in the simplified 1 + 1 -single photon picture. 

5.1 .1 . Distribution of tine liigli frequency plioton 

As an example, we consider the decay of the 3s-level in vacuum. 
If there are no perturbing particles, two photons will be released 
and the emission profile (see Fig. O is described by the two- 
photon formulae discussed in the previous sections. In the 1 -i- 
1 -single photon picture, with very high probability the electron 
after a short time (~ 1.6 x 10-^ s) decays to the 2p-state, emitting 
a photon close to the Balmer-ff frequency. Then it independently 
releases a second photon, for which the frequency distribution, 
in the rest frame of the atom, is given by the natural line profile. 
Therefore the number of photons appearing per second in the 
frequency interval v and v + dv in the vicinity of the Lyman-a 
resonance due to the transition from the 3s-state is given bjOB 



6- , '(v) dv = 

2p^ls V / 



•2p 



/ply 

^ 2p^ls 



-,V-Va 



dv. 



(16) 



where L{a, b) is defined in Eq. (fT2l i and v„ is the Lyman-a central 
frequency. Integrating over v yields / <^2p^Ys'^*-^-' ~ ^3s-^2p' 
i.e. the total rate at which electrons are added to the 2p-state 
after the transition from the initial 3s-level. Note that here we 
have assumed fap ~ ^^2p^is- 

For (fTST i one assumes that there is no coherence or corre- 
lation between the first and second photon, and consequently 
the Lyman-a line -profile in the 1 + 1 -photon picture is a pure 
Lorentzian up to very large distances from the resonance. This is 



The factor of 1 jn is required due to the normalization of L{a, b). 
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also equivalent to assuming that the transition from the 3s-state 
leads to a 'natu ral' distribution of electrons within the 2p-state 
dMihalas'igTSV One can also obtain this result using the inter- 
pretation of Weisskopf & Wigner ( 1933) for the line width. 

Looking at other initial s or d-states, the same argument as 
above can be carried out. In the more general case, one simply 
has to replace AjJ^jp with the corresponding partial spontaneous 

decay rate A^^^^^2p '■^^ 2p-state. This shows that no matter 
what is the initial level, the shape of the 1 + 1 -emission profile 
in the vicinity of the Lyman-a resonance is always a Lorentzian. 
Within the 1 + 1 -single photon picture, the same is true for the 
other possible intermediate resonances (e.g Lyman-yS, y, etc) in 
the two-photon cascades from high initial s or d-states. However, 
there in addition the partial width of the 2p-state due to the tran- 
sition to the ground level appearing in Eq. (fTSI l has to be replaced 
by the corresponding total (one-photon) width of the intermedi- 
ate p-state. Also one has to take into account the branching ratio 
for transitions leading directly to the ground state. 
With these comments one then can write 



1 + 1 -single photon formulation: 



^ly ply 
f> , (v) dv = ; L 



4^ 



, V - V„is 



dv, (17) 



where v„is is the central frequency of the corresponding Lyman- 
series transition. 



5.1 .2. Distribution of tine low-frequency plioton 



In Sect. 15.1.11 we have focused on the high-frequency photons 
released in the two-photon cascade. If we now consider the low- 
frequency photons, then the profiles of these will be given by 



iiy 



(v) dv = 



71 



/ply 
4^ 



, V - v„,, 



dv, (18) 



where v„ „p - v„ is - v„is is the transition frequency from the 
initial «/S or n,d-state to the np-state. 

Here one may ask why the width of the line is determined by 
the width of the intermediate «p-state only and not by ^n^+^J^i^ 
as usual. We simply wanted to be consistent with the approxi- 
mate treatment of the cascade and interference terms in the full 
two-photon formulation (see Sect. 12.1.31 ). for which the width 
of the initial state was neglected. As mentioned above, physi- 
ca lly this is motivated by the idea that, within the formulation 
of IWeisskopf & Wignej i[T930l) . we consider one particular ini- 
tial 'energy sub-level' and do not specify the process that pop- 
ulated it. Therefore the final profile is independent of the shape 
of the distribution of energy-sub-levels around the mean energy 
of the initial state. One can also consider this as equivalent to 
neglecting any possible reshuffling of the electron by perturbing 
particles while it is in the initial state. Furthermore, in general 
r'n > r'^ , , such that r'^ , , would not contribute much to the 
total width of the line. But most important, in our computations 
we do not approach the resonances so close that these differences 
would play any role. 



5.1.3. Total two-photon profile in the single photon picture 

With Eqn. ( fTTT i and ( fTSl l it is now possible to write the total dis- 
tribution of photons emitted in the two-photon decay of an iso- 
lated hydrogen atom in some given initial s or d-level within the 



"n,s/d- 



.isM-ZKX,pM + C^fi:"''Ml- (19) 



;i=2 



Integrating 



j^(v) over all possible frequencies yields 



i r 0'+'^ , (v)dv = 2:"'VA'^. wWch is the total 

contribution to the width of the initial level due to spontaneous 
transitions that lead to the ground state and involve exactly two 
photons (cf. Eq. Q). 

What are the main differences of the 1 4- 1 -photon profile with 
respect to the full two-photon profile, as defined by Eq. ^1 

(i) There is no non-resonant contribution, resulting from virtual 
transitions via intermediate states with energies higher than 
or equal to the initial state. 

(ii) As a consequence of (i) there is no resonance/non-resonance 
interference term. 

(iii) In contrast to Eq. (I14ab . there is no interference among the 
resonances. 

(iv) As a consequence of (ii) and (iii) each resonance has the 
shape of a Lorentzian up to very large distances from their 
line centers. 

(v) Usually one does not restrict the range of integration to the 
interval < v < v„is. 

It is also important that interpreting each resonance appearing in 
^]vsid^\S^^ separately, it is possible to uniquely define the rates 
at which electrons flow in and out of a particular intermediate p- 
state. In astrophysical computations this is the usual procedure 
for solving the radiative transfer problem for each transition sep- 
arately. Within the full two-photon formulation this is only pos- 
sible very close to the centers of the resonances (where the con- 
tribution from the other terms is negligible), but in the wings, 
photons from different intermediate transitions contribute non- 
trivially and make this separation difficult. Also this overlap of 
the resonances is taken into account in Eq. ( fT9] l, but it is usually 
neglected in astrophysical computations. 

As an example, we illustrate the differences in the two- 
photon emission profiles for the initial 4s and 4d states in Fig. 
|5] One can see that in the distant wings of all the resonances 
the differences of the profiles are rather big. This is mainly due 
to the non-resonant term and its interference with the cascade 
contribution, but also the resonance/resonance interference plays 
some role. Below we now focus our analysis on the deviations 
of the two-photon profile from the pure Lorentzian close to the 
Lyman-ff resonance. These differences are the main reason for 
the corrections to the hydrogen recombination history. 



5.2. Two-photon profiles close to the Lyman-a resonance 

In low-density plasmas like the expanding Universe during 
cosmological hydrogen recombination, har dly any perturb- 
ing particle is within the Weisskopf-radius (I Weisskopj 119321; 
ISobelmanetanil995h . so that the coherence in two-photon and 
possibly multi-photon transitions is maintained at least for the 
lower shells. In astrophysical computations the frequency dis- 
tribution of photons released in the Lyman-a transition due to 
electrons reaching the 2p-state from higher levels is usually 
described by a pure Lorentzian. Within the interpretation of 
IWeisskopf & Wigneil(ll930l) . this means that the electron is com- 
pletely reshuffled among all the possible 2p energy sub-levels. 

In calculations of the cosmological hydrogen recombination 
problem, we are now interested in the deviations of the full 
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Figures. Comparison of the two-photon emission profiles for the 4s and 4d states. We show, <('^^/^_^i^, following from the full 2j 



treatment according to Eq. (|7]i, and, il>\^!j_^i^, using the 1 + 1-single photon description, as given by Eq. (T9[ . The first 200 terms 
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Figure 6. Normalized 1 + 1 -two-photon profile 02p-^is^ ^' 'US' close to the Lyman-or frequency in comparison with the 

full two-photon-profiles, ^j^s^j^js according to Eq. (|7l), for several initial s and d-states. We divided the two-photon spectra by their 
partial one-photon transition rate to the 2p-state (r„s/d^2p) and transformed to the variable ja - vlVa- For the given curves we 
included the first 500 terms above «,. The dashed curve corresponds to the usual Lorentzian profile of the Lyman-a transition. 



two-photon profile from the normal Lorentzian shape. Here 
one should mention that in general the deviations of the 2y- 
emission profile, using the full two-photon treatment as de- 
scribed in Sect. |2l from the one in the 1 4- 1-single photon de- 
scription (Eq. ( fT9t ) should also be considered. However, in the 
red wing of the Lyman-a resonance, one can write 



^n;s/d^ls 



Z 

n= 

A 



n=2 



,17 



/((S/d— >np np— >ls 



4;r2 



[v - V„,„p] 



47r2 



[v- V„ls] 



,ir 



H,s/d— >2p 2p— >ls 

47r2 [v - V2pis]^ 



(20) 



where the last step is possible, since the distant wings of all 
the other resonances do not contribute significantly when one is 
close enough to the center of the Lyman-o' transition. The more 
one approaches the Lyman-ff resonance, the better this approx- 
imation becomes. For the estimates carried out below, this ap- 
proximation is satisfactory. 

To understand the deviations of the two-photon emission 
profiles close to the Lyman-a resonance, we now directly com- 
pare according to Eq. (|7]i with 02p^is^'^ given by 
Eq. ([TtI i. For convenience we choose ja - v/v„ as the common 
frequency variable. Then the full two-photon profile in this new 
coordinate is given by (p^iya) - Va(l\^/^^i^{y)/viu. The axis of 
symmetry is then aty,, - |[«^ - instead of y = 1/2. Since 
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0.99 



0.99 



= V / V 



= V / V 



*nI/d^i,s(^'>-)/'^"is - 'f'2j^ls^^(ya)]/(f>2j^u^^(ya), of the curvcs given in Fig.|6]with respect to nor- 



Figure 7. Relative difference, 
malized 1 + 1 -two-photon profile. Very close to ya ~ I, the curves should be considered as extrapolated estimates. 



in the vicinity of any particular resonance all the two-photon pro- 
files scale like r„s/d^„'p, focusing on the Lyman-a transition, we 
also re-normahzed by r„s/d^2p- 



In Fig.|6]we give the normalized 1 + 1 -two-photon profile in 
the vicinity of the Lyman-ff transition in comparison with the re- 
normalized two-photon-profiles for several initial s and d-states. 
One can see that at large distances the two-photon profiles in 
the full 2y-treatment deviate a lot from the Lorentzian shape. 
For both the initial s and d-states, the very distant red wing is 
several times above the Lorentzian. Within the frequency range 
0.9 < ya ^ 1-1 for the s-states, the red wing lies below, the 
blue wing above the Lorentzian, whereas the opposite is true for 
the d-states. In particular for the d-states, the red wing is always 
above the Lorentzian, and unlike the s-states in the considered 
frequency range there is no additional zero below the Lyman- 
a resonance. In Fig. |6]one can also see that for the chosen set 
of coordinates, the variation in the profiles is rather small in the 
case of initial s-states, and the modifications become negligible 
even for initial d-states above « ~ 10. 



In Fig.|7]we show the relative difference of the curves given 
in Fig. |6]with respect to the Lorentzian of the Lyman-a reso- 
nance. The wing redward of the Lyman-o- frequency lies below 
the Lorentzian for initial s-states, exceeding the level of ~ 10% 
at more than ~ 1 .6 x 10'* natural width from the center. For initial 
d-states, in all shown cases the wing redward of the Lyman-a 
frequency lies above the Lorentzian. The relative correction to 
the Lorentzian scales roughly linearly with Av = v - Vq. in this 
regime. Therefore the net change in the rate of photon produc- 
tion in the red wing of the Lyman-a transition at frequencies in 
the range vi < v < V2 depends logarithmically on the ratio of vi 
and V2: ANy oc ' ^ x Avdv oc log[v2/vi]. Here we used 

the wing approximation of the Lorentzian L oc l/[v - v„f. This 
estimate shows that the value of the effective two-photon decay 
rate does not depend very strongly on V2 (see Sect. l5.3l l. 



5.3. Cosmological hydrogen recombination 

5.3.1 . Escape of photons in the red wing of the Lyman-a 
resonance 

In the context of cosmological hydrogen recombination, the 
escape of photons in the red wing of the Lyman-a res- 
onance, which is one of the major channels to reach the 
ground state of hydrogen, pl ays a key role in controlling the 



dynamics of recombination fVarshalovich & Svu nvaevI 1 19681 : 
Grachev & Dubrovich 1991; Rybicki & dell 'Antoniol 1 19941: 
iChluba & Sunyaevll2007aft . At large distances, say at frequen- 
cies below Vc redward of the Lyman-a central frequency, Va, the 
probability of absorbing a photon to the continuum, thereby cre- 
ating a free electron, becomes very low. Photons released below 
Vc directly escape further interaction with the neutral hydrogen 
atoms and lead to the settling of an electron in the Is-state. On 
the other hand, all photons emitted at frequencies v > Vc will 
have a very high probability of being absorbed in the contin- 
uum or undergoing transitions to higher levels, possibly after 
many interactions with neutral hydrogen atoms or when redshift- 
ing into the domain of the Lyman-a resonance from frequencies 

V > Va. 

Determining the exact value of Vc during the epoch of cos- 
mological hydrogen recombination requires a full treatment of 
the radiativ e transfer in the Lyman-a r esonance. Our computa- 
tions show dChluba & Sunyaevll2007ah that Vc depends on red- 
shift and should typically lie within 100 to 1000 Doppler width 
below the Lyman-a frequency. At redshift z, one Doppler width 
coiTesponds to 



Avn ^ 58.0 



92.5 



(1+z) 



1100 

(1 +z) 



1/2 



GHz ^ 2.35 X 10" 



(1+z) 



1100 



1/2 



1100 



1/2 



2p^ls' 



(21) 



hence |vc - Vq,|/v„ ~ 10 ^ - 10 ^, or in terms of the 2p-line width 



■yj/r. 



ir 

2p^ls 



10"* - 10^ 



In computations of the hydrogen recombination history, it is 
therefore important to know how many photons reach the very 
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3d -> Is 
4d -> Is 
5d "> Is 
6d -> Is 
8d -> Is 
10d--> Is 



r = [v-v ] / r. 



r = [V-V ] / 

r a-" 2p 



Figures. Rate of photon production at frequencies below xp according to Eq. ( |22] | for several initial s and d-states. To convert to 
the variable [v - Vq,]/v„, one should multiply xp by ~ 2.54 x 10"^. 





5000 



r = [V-V ] / r. 



otil 

-10' 



r = [V-V ] / r. 



Figure 9. Effective change in the rate of photon production (real profile minus Lorentzian) at frequencies below xp, according to 
Eq. (l24l i. To convert to the variable [v - v„]/vq,, one should multiply xp by ~ 2.54 x 10"^. 



distant red wing of the Lyman-ff resonance directly. If we want 
to estimate this effect, we need to compute the difference in the 
number of photons, that are directly escaping in the distant wing 
by comparing the emission profiles in the full treatment of two- 
photon processes with the one in 1 -H 1 -single photon picture. 
This will show the relevance of this process. 

If we consider those photons emitted in the red wing of the 
Lyman-a resonance because of two-photon transitions from up- 
per s or d-states, then when introducing the dimensionless fre- 
quency variable xp(v) = [v - VQ,]/r2p, the results discussed in 
Sect. 152] suggest the following: 

(i) Because of two-photon processes, more photons will escape 
in the very distant red wing of the Lyman-a resonance (be- 
low xp 5 X 10^) than in the 1 + 1-single photon treatment 

for the direct cascade emission. 

(ii) For initial s-states, fewer photons are emitted in the range 
-5 X 10^ < Xp < -10"* than in the 1 + 1-single photon picture. 



(iii) For initial d-states, more photons are emitted in the range 
-5x10^ < Xp < -10"* than in the l-H 1-single photon picture. 

Because of (i) and (iii) hydrogen recombination should occur 
slightly faster, while (ii) may make it a bit slower. Since the sta- 
tistical weight of d-states is 5 times higher than for s-levels, one 
expects that recombination will in total be slightly faster than in 
the standard treatment because of two-photon processes. 

In addition to the direct escape of photons in the distant red 
wing of the Lyman-or transition, also significant differences close 
to the line center arise (see Fig. |7]i. Understanding how these 
changes affect the effective escape of photon from the line center 
requires a more rigorous treatment of the radiative transfer prob- 
lem in the line. Also the feedback of photons emitted in the blue 
wing of the Lyman-ff transition and in particular those coming 
from the other Lyman-series transitions, should be slightly mod- 
ified when taking the full two-photon process into account. Both 
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aspects are beyond the scope of this paper and will be addressed 
in a future work. 

We can now estimate the effect of the changes in the effective 
escape of photons in the distant red wing of the Lyman-o' transi- 
tion. For this only the photons between the innermost resonances 
in the two-photon emission spectrum are contributing (e.g. pho- 
tons between the Balmer-y and Lyman-c transition for the 5s and 
5d-two-photon decay, see Fig.[3]i. This is because we only want 
to count photons up to v < and correspondingly v' > v„is - Vc- 
Because of the symmetry of the full two-photon profile, it is 
therefore sufficient to integrate <;*„s/d^is^^) from y - 1/2 up to 
Jc = Vc/v/is: 



^Jy)dy. (22) 



This integral yields the total number of photons that directly es- 
cape per second in the distant red wing of the Lyman-o- line. It 
should be compared with the value computed using the standard 
1 + 1 -single photon profile. Since we only consider cases very 
far in the red wing of the Lyman-a transition, the integral over 
the Lorentzian resulting in the 1 + 1 approach can be written as 



.«s/d,l+ly. s 



Jo 



rflS/d-l+lV/ N J 



^ns/d^2p 1 + Xr,c£ 



47r2 



(23) 



where we introduced e = r2p^is/v„ ~ 2.540 x 10"^ and used 
the variable xr(v) = [v - v„]/r2p. The effective difference in the 
photon production rate, or equivalently the photon escape rate in 
the distant wings, is then given by 



*2p^ls 



(24) 



for a fixed frequency v^. Although in general Vc is a function 
of time, below we assume that it is constant. A more rigorous 
treatment will be presented in some future work. 

5.3.2. Approximate inclusion into tine multi-level code 

In our formulation, AAns/ri^is(Vi.) plays the role of the pure two- 
pho ton rate coeffic i ents u sed in iDubrovich & GrachevI ( |2005|) 
and I Wong & Scot J (l2007h . If we want to estimate the possible 
impact of our results on the hydrogen recombination history, we 
have to take the additional net escape of photons into account. 
This can be accomplished by adding 



gw-s/d 
gls 



(25) 



to the rate equation of the Is-state and subtracting it from the 
corresponding rate equation of the «s and «d-levels. Here A^is, 
N„s, N„d are the number density of hydrogen atom in the Is, 
«s, and «d-states, respectively. Furthermore, Ty = ro(l + z), is 
the temperatu re of the ambient black body radiation field, with 
To = 2.725 K dFixsen & Mathed2002h . The factors g„s = gu and 
grid = 5 gu are due to the statistical weights of the s and d-states. 
In ( |25] ) we have neglected any possible deviation in the radiation 
field from a blackbody and also omitted stimulated two-photon 
emission. Both processes should only lead to higher order cor- 
rections. Moreover we have added an inverse term, assuming de- 
tailed balance. This term is not important during the main epoch 
of hydrogen recombination (z < 1600) and was only included 
to maintain full thermodynamic equilibrium at high redshifts. A 
self-consistent derivation is beyond the scope of this paper. 



5.3.3. Results for the photon production in the distant wings 

In Fig. [8] we give the rate of photon production at frequencies 
below xr within the full two-photon treatment, i.e. according to 
Eq. ( l22b . for several initial s and d-states. For the d-states the 
photon production is ~ 10 times faster than for the correspond- 
ing s-state. In the case of initial s-states, the plateau of ^JJs_^is(vc) 

close to Xr 4 x 10^ is caused by the zero in the central region 

of the two-photon emission spectra (e.g. see Fig. |6l). As men- 
tioned in Sect. [3] this zero is absent in the two-photon spectra 
of initial d-states, and consequently no such plateau appears for 
A^^^j^(Vc). In both cases the rate of photon production decreases 
when increasing n. Looking at Fig. H] just from the non-resonant 
term one would expect the opposite behavior However, due to 
destructive interference this does not happen. 

In Fig. |9] the net change in the rate of photon production 
at frequencies below xr is shown. The photon production due 
to the two-photon decay of initial s-states, at relevant distances 
from the Lyman-a center (jcp > -10^), is actually slower than 
in the 1 + 1 -single photon picture. This suggests that due to the 
full treatment of the two-photon process for the s-states alone, 
cosmological hydrogen recombination is expected to be slower 
than in the standard c omputations. This contrasts to the work of 
IWong & Scotll (|2007|) . where both the s and d-state two-photon 
process leads to an increase in the rate of recombination. 

On the other hand, for the d-states the effective photon es- 
cape rate is higher than in the 1 -i- 1 -single photon picture, hence 
one expects an increase in the rate of recombination. Since the 
statistical weights of the d-states are 5 times larger than the s- 
states, and also the effective increase in the wing photon produc- 
tion rate is roughly additional 5 times higher (cf. Fig. |9]l, one 
still expects that, even when including the combined effect of 
the s and d-state, two-photon process, cosmological hydrogen 
recombination in total will proceed faster than in the standard 
treatment. 

We would like to mention that using the analytic approxima- 
tions given in AppendixiBlfor the non-resonant term in connec- 
tion with the formulae in Sect. 12. 1.3] we were able to reproduce 
the rates presented in this section. 



5.3.4. Comparing with earlier works 

In the works of lDubrovich & GrachevI (l2005h and lWong & Scottl 
(|2007|) . only the combined effect of the two-photon process for 
the ns and nd-states on the hydrogen recombination history was 
discussed. To compare our results for the effective photon pro- 
duction rates with their values, we also write the combined ef- 
fective decay rate 



6AA„ 



is(vc) = AA„s^is(vc) + 5 AA„ 



s(Vc), 



(26) 



where we implicitly assumed that the ns and nd-states are in full 
statistical equilibrium with each other {N„d = 5 A^„s). At the rel- 
evant redshif t, the deviations from full statistical equiUbrium are 
rather small dRubino-Martm et al.ll2006t IChluba et all2007h . so 
that this approximation is possible (see Sect. l5.3.5T l. 

In Fig. [TOlwe present the results for AA„s+„d^is(vc) for sev- 
eral shells. If we consider the effective rate for the 3s and 3d- 
levels then, even for very conservative values of Vc, say 1000 
Doppler width or ~ 10^ natural width below the Lyman-o- reso- 
nance, we_obtainAA^s+M^^>is^ w hereas from the formu- 
lae in IDubrovich & GrachevI (l2005h and IWong & Scotl (l2007l) 
one can find AA?"^,. , ~ 22 s"' and AA*^, , , ~1.5s-\re- 

3s+3d^ls 3s+3d^ls ' 

sp ectively. Our value is only ~ 3.3 times lower than the one 
of lDubrovich & GrachevI (l2005l) but ~ 4.5 times higher than in 
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3s+3d-> Is 
4.s+4d-> Is 
5s+5d -> Is 
6s+6d — > Is 
8s+8d-> Is 
l()s+l()d-> Is 



-10" -10- -10' -^"^ 

jt = [v-v ] / 

Figure 10. Combined effective two-photon photon production 
rate, AA„s+„d^is(i'c), computed according to Eq. ( |26] |. To con- 
vert to the variable [v - Va]/va, one should multiply xp by 
~ 2.54 X 10-^ 

Table 2. Effective difference in the photon production rate in the 
distant wings using Eq. (l24l i. 



n 


^r,c =-10= 


xr.c = -5 X 10* 




3 


-0.307/8.133 


-0.951 /9.882 


-2.565/ 13.883 


4 


-0.210/4.619 


-0.461 /5.648 


-1.082/8.004 


5 


-0.122/2.523 


-0.245 / 3.082 


-0.551 /4.356 


6 


-0.074/1.497 


-0.144/ 1.8272 


-0.317/2.578 


7 


-0.048 / 0.954 


-0.091/1.164 


-0.199/ 1.641 


8 


-0.033 / 0.643 


-0.061 /0.785 


-0.133/1.106 


9 


-0.023 / 0.454 


-0.043 / 0.553 


-0.093 / 0.779 


10 


-0.017/0.331 


-0.032 / 0.404 


-0.068 / 0.569 



IWong & Scot3 (|2007|) . We would argue that for the third shell 
even values up to 10 s"' still are reasonable, in particular at 
very low (z < 1000) and high (z > 1300) redshifts, where 
the probability of absorption decreases. In Table |2] we give a 
few values of AA„s/d^is(Vc) for different frequencies v^. Given 
are the values of AA„s/d^is in 1/sec for different frequencies 
■'^r.c - [vc - Vff]/r2p^is- In each column the first value is for 
the s-levels, the second for the d-states. 

Figure [TOl also shows t hat, in contrast to the w orks of 
iDubrovich & Grachev! ( l2005h and lWong & ScottI ( |2007|) . the net 
photon escape rate due to the combined effect of the s and d- 
state, two-photon process decreases with increasing n. This im- 
plies that the relevance of the two-photon emission from higher 
shells is significantly less than in their computations, because the 
sharp drop in the populations of levels with n will no longer be 
partially canceled by the assumed linear increase in the effective 
two-photon-decay rate. 

5.3.5. Differences in tine free electron fraction 

We modified our multi-level hydrogen code (for more details 
see iRubino-Martm et all 120061: Ichluba et alj|2007l) to take into 
account the additional escape of photons in the distant wings 
of the Lyman-a resonance due to the two-photon process us- 
ing Eq. ( |25] ). For the hydrogen atom we typically included the 
first 30 shells in our computations, following the evolution of 



the populations for each angular-momentum substate separately. 
We also performed computations with more shells, but this did 
not alter the results significantly with our approach. 

The additional two-photon process was included for s and 
d-states with n < n2y, where the parameter n2y gives the high- 
est shell for which the additional two-photon decay was taken 
into account. We only used «2y ^ 10, but because of the strong 
decrease of |AA„s/d^is(Vc)| with n (see Fig. |9]l and the drop in 
the populations of higher shells, we do not expect any signifi- 
cant differences when going beyond this. For simplicity we also 
assumed that the value of Vc is constant with time. This makes 
our estimates more conservative, since both at very low and very 
high redshifts, Vc should be closer to v„ and therefore may in- 
crease the impact of the two-photon process on the recombi- 
nation history. We performed computations with three different 
values of v^. The effective rates for these cases are summarized 
in Table|2l We consider the case with jcp.c - -10^ as pessimistic, 
whereas the case xp.c = -10^* may be optimistic. 

We also ran computation s usi ng the formula e accor ding to 
IDubrovich & GrachevI (lIOOSi) andlWong & ScottI (l2007l) . In the 
paper of IDubrovich & GrachevI (l2005h . the s and d-rates were 
not given separately, but assuming AA^^^ = a 4 
plicity, one finds 



, = AA" , , for sim- 



AA'^-°P = AA^l-""? = 8.2293 s" 



xSocin), 



(27) 



where 5dg(«) 



-41 



Since the deviations from 



full statistical equilibrium are rather small dRubino-Martm et"aLl 
12 061 ; IChlub"a et al.ll2007h . this assumption should not be very 
crit ical and, in any c ase, is only meant for comparison. 

IWong & ScottI ( |2007|) explicitly give the rates for the 
3s and 3d-states and then assume the same «-scaling as 



IDubrovich & GrachevI (l2005l) . This yields 

AA„'^:^^ = 8.2197 s-ix5ws(«) 
^'d^Tfs^ 0.13171 s-'x5ws(«), 

with5ws(«) = 5dg(«)/5dg(3) = f [f^] 



(28a) 
(28b) 

Comparing 

wi th Eq. (p7|). one can see that the d ifference in the appro ach 
of IDubrovich & GrachevI (l2005l) and IWong & Scotl dioOTl) is 
mainly because they used a much lower rate for the d-states (a 
factor ~ 170!). The assumed rate for th e s-states is only ~ 2.7 
times lower than in the computations of IDubrovich & GrachevI 

mm- 

In Fig.[TT|we present the relative change in the free electron 
fraction when only including the additional two-photon process 
for 3s and 3d-states. For co mparison we show the results ob- 
tained using the decay rates of lDubrovich & GrachevI (l2005b and 
IWong & Scoti (120071) . One can clearly see that the dependence 
on the adopted value of Vc is not very strong. For our optimistic 
value of Vc, close to the maximum th e effect is roughly 2 times 
smaller than for the values of Dubrov ich & Grach ev ( 2005i), and 
even in our pessimistic model, it is stil l more than ~ 4 time s 
greater than within the framework of IWong & Scota (I2007h . 
C omparing the curves, which w e o btained within the app roach 
of lDubrovich & GrachevI (120 05') and'Wong & ScottI ( 120071) . with 
those in Fig. 3 of lWong & Sc ott ( 2007) one can see that our re- 
sults for the changes in the electron fraction are slightly smaller 
We checked that this is not due to our detailed treatment of the 
angular-momentum substates. This is expected since the devi- 
ations from full statistical equilibrit im at the relevant redshifts 
are too small to hav e any effect here dRubino-Martm et a ni2006l : 
IChluba et al.ll2007l) . Also we computed the same correction us- 
ing 50 shells, but found no significant increase. 
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Figure 11. Relative change in the free electron fraction. Here 
we only included the additional two-photon process for the 3s 
and 3d-states. The computations were performed for a 30-shell 
hydrogen atom. The effective two-photon rates for three differ- 
ent values of Vc according to Table |2] were used. For compar- 
ison, we show the results th at were obtained using the 3s and 
3d two-photon decay rates of iDubrovich & Grachev! (l2005l) and 
IWong&Scot1(l2007h . 



Figure 13. Relative change in the free electron fraction when 
taking the additional two-photon emission for up to 10 shells 
into account. The computations were performed including 30 
shells, for three different values of Vc. For comparison, th e resul t 
that was obtained within the approach of IWong & ScottI (|2007|) 
is shown, but for «2y = 40 and using a 50-shell model for the 
hydrogen atom. 




Figure 12. Relative change in the free electron fraction for dif- 
ferent values of «2y ■ The computations were performed including 
30 shells. 



In Fig. [12] we illustrate the impact of the two-photon pro- 
cess from higher shells. With our estimat es of the effective two- 
photon dec ay rates, like in the studies of IDubrovich & GrachevI 
(l2005h and lWonp & Scott (2007), the effect increases with niy 
However, the strong decrease in the effective rates within our 
computations (see Table |2]i implies that the result practically 
does not change when including the additional two-photon ef- 
fect for more than 5 shells. T his s trongly contra sts the works of 
IDubrovich & GrachevI ( l2005h and lWong &"Scott (2007), where 
the total change in the free electron fraction radically depends 
on the chosen value of «2y (even up to «2r = 40 was consid- 



ered). As mentioned above, in these computations the increase 
in the two-photon decay rates with n (cf. Eqn. dZTI i and ( l28T l) 
partially cancels the decrease in the population of the higher lev- 
els, and therefore enhances the impact of their contribution as 
compared to the lower shells. For example, at z ~ 1200 (i.e 
close to the maximum of the changes in A^e) the populations 
of the excited states are still nea rly in Saha-equilibrium with 
the continuum dChluba et alJl2007l) . Therefore the population of 
the fourth shell is roughly a factor of ^ '"'41/^^)' ~ 0.19 smaller 
than in the third shell. Also the eff ective 2y-ra te decreases by 
~ 1 .8, whereas in the pic ture of D ubrovich & Grachev ( 200$) 
and lWong^ & ScottI (l2007h it would have increased ~ 1.9 times. 
From Fig. [I2|it is also clear that the strongest effect for our es- 
timates of the effective decay rates comes from the 3s and 3d- 
levels a lone. This again is in strong opposition to the computa - 
tions of IDubrovich & Grachev (2005) and Wong & Scott! (|2007|) 
where more than ~ 75% of the correction is due to the combined 
effect of higher shells. 

In Fig.[T3]we give our final estimates for the possible changes 
in the recombination history. In our optimistic model the change 
is AXe/Xe ~ -0.53% at redshift z ~ 1150, and it drops to 
AXs/Xg ~ -0.32% for the pessimistic case. Including more 
shells in the model for the hydrogen atom did not change these 
results. For comparison we also computed the changes in the 
ioniza tion history by applying the formulae of IWong & ScottI 
(I2OO7I) . but using n2y = 40 and 50 shells for the model of the 
hydrogen atom. Alt hough our discussion has shown that the 
values computed by ICresser e"ta D (119861) for 3s and 3d-states 
are likely not related to the cosmological hydrogen recombi- 
nation problem and that extrapolating those values to higher 
shells is rather rough, our final results are nume rically compat- 
ible w ith those obtained using the approach of IWong & ScottI 
(120071) . However, within our approach the changes in the ion- 
ization histo ry close to the maximum of the Thomson visibil- 
ity function dSunvaev & Zeldovich| 19701) are larger than in the 
computations of IWong & ScottI (I2007h . Therefore the changes 
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in the cosmic microwave background temperature and polariza- 
tion power spectra are also expected to be a bit larger Indeed 
it seems that the corrections to the ionization history due to the 
two-photon decay from higher shells does not reach the percent 
level, and that the impact of this process was overestimated by 
iDubrovich & GrachevI (l2005h . 

5.3.6. Additional remarl<s 

Here we have only investigated the bound-bound two-photon 
transitions directly leading to the ground state. Equations ([T]i and 
^ are also applicable when the final state is any s-level. We 
also checked the rate for the two-photon transition 3s — > 2s and 
3d — > 2s and, as expected, found very low values (~ 0.0885 s"' 
for the 3s and 0.0278 s ' in the case of 3d). In addition, be- 
cause all dipole transitions to the second shell are optically thin 
in the recombination problem, these corrections should never 
be important within this context. Similarly, the 2p — * Is two- 
ph oton transition due to it s low probability (~ few x 10"* s"', 
see iLabzowsky et aPfeOOSh can be completely ignored. 

One may in addition consider the problem of two-photon 
transitions starting from the continuum, e.g. the recombination 
of electrons to the 2p-state and subsequent release of a Lyman- 
a photon. Here deviations of the line profile from the normal 
Lorentzian shape can also be expected and may lead to an in- 
crease in the effective Lyman-ff escape rate. However, since the 
supply of photons to the 2p-state by transitions from higher 
shells is several times faster, the total impact of this effect is very 
likely less than the one from the 2y-transitions already discussed 
here. 

As mentioned above, under physical conditions like those in 
our Universe during the epoch of hydrogen recombination, the 
coherence of two and possibly multi-photon processes is main- 
tained. Consequently one should investigate how strong the devi- 
ations of the corresponding emission profiles from a Lorentzian 
are when more than two photons are involved. This requires a 
QED multi-photon treatment, which is beyond the scope of this 
paper But as we have seen above, adding the two-photon pro- 
cess for the 4s- and 4d-states (due to the drop in the population 
of these levels and decrease in their effective two-photon decay 
rate) has aff'ected the recombination history at a level of ~ 0.02% 
in addition to the 3s and 3d (see Fig.[T2]). For the 3y-decay of the 
4f-level, one expects that the relative correction will be close to 
the one from to the 3d two-photon decay. This is because the 
largest term in this 3y-description should involve at least one 
nearly resonant transition, since the other contributions should 
be suppressed in addition. For photons appearing close to the 
Lyman-a line, a nearly resonant 4f — > 3d transition, followed by 
a quasi 2'y-decay of the 3d-state, is most likely. Therefore one 
has 

where we have usecO A4^_^3j ~ ^if_^3d and Ag^^^p ~ ^3d^2p- 

With the value in Table |2l one finds AA'^f^^^ ~ 2.1s-'. Due 
to the statistical weight of the 4f-state, one then should obtain 
another ~ 7x2.1s-' ~ 14.7s-' addition to the overall decay 
rate of the fourth level. This is roughly a ~ 50% correction to 
the total two-photon correction of the fourth shell. Therefore, 
one could expect another ~ 0.01% correction to the free elec- 
tron faction when including the three-photon decay from the 4f. 

' ' In both cases the branching ratios are very close to unity. 



Furthermore, for some high level (likely above n ~ 20-50), even 
the residual disturbances from the small amount of perturbers 
that are present in the Universe will destroy the coherence of the 
possible multi-photon emission coming from very high initial 
levels. 

One could also consider the three-photon decay of the 2p- 
state. Here, just as in the 2s-two-photon decay, no intermediate 
resonances are involved; and due to momentum conservation, 
this process is allowed. However, simple estimates show that this 
process has a rate lower than ~ 10"^ s ' and hence negligible at 
the 0.1%-level. 

6. Conclusions 

We have studied in detail the emission of photons due to two- 
photon transitions from high s and d-states to the ground level. 
Up to n - 20 we found simple analytic fitting formulae to rep- 
resent the full two-photon emission profile with very high accu- 
racy. We have discussed the deviations in the two-photon emis- 
sion profiles from the natural Lorentzian shape and investigated 
the importance of the non-resonant, cascade, and interference 
term separately. 

Applying our results to the cosmological hydrogen recom- 
bination shows that the corrections to the ionization history due 
to the additional two-photon process from higher shells likely 
do not reach the percent level. For conservative assumptions we 

find a correction AX^/Xe 0.4% at redshift z ~ 1160. This is 

numerically similar to the result of I Wong & Scotll ( |2007|) : how- 
ever, the physics leading to this conclusion is rather different. 
In particular we find that the two-photon process for initial s- 
states actually slows the recombination process down. In addi- 
tion, the effective two-photon rates connecting the high s and 
d-level directly to the 1 s-level decrease with principle quantum 
number n. Both aspects contrast to the rate estimates used in 
the stu dies by Dubrovich & Grachev (2005) and Wong & Scotll 
(I2007h . Here it is very important that the destructive interference 
between the cascade and non-resonant term cancels a large part 
of the additional non-resonant emission in the distant red wings 
of the Lyman-ff transition. Furthermore, in our computations the 
main correction to the ionization history stems from the 3s and 
3d-sta tes, while in the computat ions of IDubrovich & GrachevI 
(l2005h and lWong& Scotll (l2007h more than ~ 75% of the cor- 
rection is due to the combined effect of higher shells. 
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Appendix A: Radial integrals 

A.1. Bound-bound radial integrals 

For the the required bound-bound radial integrals up to n, = 5, one has 

(n - l)"-5/2 



{R. 
{R. 
{R. 
{Rs. 

{Rm 

{R4d 
(«5d 



r\R„p 
r\R„p 
r\R„f 
r\R„p 
r\R„p 

r\R„p 
r\R„p 
r\R„p 



:2'n 



(n + l)"+5/2 
:2'*V2«'/2(„2_i)l/2 



(n + 2)''+3 



:2*33V3„'/2(„2_i)l/2[7„2_27] 
2" 



in - 3)"-* 
{n + 3)"+'* 

(n - 4)"-^ 



n"2(;!- - 1)"2[23(/ - 288(1^ + 768] ■ 
3 ' ' ' ' (,i + 4)"+5 



2'»5* V5 



X [91«'^ - 2545n* + 20625n^ - 46875] 

1!^1V2V3^„/, ./,(«-3r:j 

VS (« + 3)"+-* 

— - i)i/2[7„2 - 48] (ilniE! 

3 V5 (n + 4)"+5 



(n - 5f 



(n + 5)« 



3 V7 



X [29n'* - 590^2 + 2625] ■ " 



(n + 5)"+6 ■ 



(A.la) 
(A.lb) 
(A.lc) 
(A.ld) 

(A.le) 
(A. If) 
(A.lg) 

(A.lh) 



In addition one needs (R„j-t \ r \R„i) = - - P. 



A.2. Bound-free radial integrals 

For the necessary bound-free radial integrals up to n, = 5, using the defin ition 
of the radial functions for the continuum states (e.g. see §36 lLandaij|| 19771) . one 
obtains 



{Rl,\r\R,i) = 2* 



-2 arctan(.i:)/.v 



{R2.\r\R,l) = 2^^^2 



(1 + .y2)5/2 Vl _ g-2^/x 
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The value of x ranges from to 00. 



Appendix B: Non-linear fitting functions 
non-resonant emission spectra 

The Tables IBT] and lR2l contain the non-linear fitting coefficients for the non- 
resonant emission spectra. The non-resonant emission spectra ai'e then given by 
C/.^lstv) = o"0')-, with cr(y) = w^'-[ao + bo v/(l h-"'* + 62 n'"^ +/:)3 w'/6 + 
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b4 vv''^ + hs + h(, w''' + h^ w''^ + hg w' + »■-)] and w = v(l - y). In this 
definition one has M^r = <T(y)/ -\]G„.i. vy^''. The first 200 terms in the infinite 
sum were taken into account. Within the assumptions, the accuracy of these ap- 
proximations should be better than ~ 0.1%. Note that oq and bo have dimension 
sec"''^, and 0"'(I/2) has dimension sec"'. 
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Table B.l. Non-linear fitting coefficients for the non-resonant /IS Is emission spectra within the frequency range 10 0.999 
for « < 20. 





2s 


3s 


4s 


5s 


6s 




-1.3974426528 x 10' 


-2.0848318929 x 10' 


-2.5412664413 x 10' 


-2.8856705303 x 10' 


-3.1567507558 x 10' 




2.5591291935 x lO"' 


5.5198239644 x 10"' 


8.0397418838 X 10"' 


4.8110524091 x 10"' 


2.5266888083 x 10"' 


bo 


6.6487585307 x 10" 


4.6179953215 x 10^ 


7.6795669026 x 10^ 


8.3598144011 x 10- 


5.6902488765 x 10^ 


bi 


-5.1211954644 x 10' 


-9.8891531378 x 10' 


-4.2588083735 x 10^ 


5.1404483119 x 10' 


4.5519082179 x 10' 


b2 


3.4686628383 x 10' 


1.0769127960 X 10^ 


7.0510361599 X 10^ 


-7.2853184204 x 10' 


-6.7551149725 x 10' 


bi 


3.8515432953 x 10' 


6.1301937755 x 10' 


-3.6120946105 x 10' 


-4.9378611706 x 10' 


-3.8327617310 x 10' 


b4 


6.7503966084 x 10" 


-4.1466090643 x 10' 


-4.4049919721 x 10^ 


4.7383284338 x 10' 


4.4401360958 x 10' 




-3 1431679488 x 10' 


-8 401 675001 5 x10' 


7 465737Q140 x 10^ 


Q Q1Q1787175 x 10' 


7 41Q5708676 x 10' 


bi 


-1.3697636757 x 10' 


7.1191755771 x 10' 


-6.5900037717 x 10' 


-1.0379836854 x 10^ 


-7.9312304169 x 10' 


h 


2.4016413274 x 10' 


-1.8695504961 x 10' 


1.8226732411 x 10' 


2.9602921680 x 10' 


2.1866188936 x 10' 


b. 


-9.5149995633 x 10" 


2.2794449683 x 10" 


-2.8692508508 x 10" 


-2.8636914451 x 10" 


-1.9806370670 X 10" 


b9 


2.6555412799 x 10" 


-4.3875761781 x 10"' 


4.4042663362 x 10^' 


3.8579054257 x 10^' 


2.6143991878 x 10"' 




2.1303295046 x 10' 


2.5955835234 x 10' 


3.0008781647 x 10' 


3.5729374807 x 10' 


4.2778194752 x 10' 




7s 


8s 


9s 


10s 


lis 


ao 


-3.3720076517 x 10' 


-3.5420358274 x 10' 


-3.6749872667 x 10' 


-3.7779791307 x 10' 


-3.8573952826 x 10' 


P 


8.5151589950 X 10"^ 


-4.2719713121 x 10-^ 


-1.3118006971 X 10-' 


-2.1932912263 x 10"' 


-2.9226372396 x 10"' 


bo 


3.7656038819 X 10^ 


2.7408956380 x 10- 


2.8136510570 x 10- 


1.5264948382 x 10^ 


8.6437136666 x 10' 


bi 


2.7980511952 x 10' 


7.4703207553 x 10" 


-1.8981579892 x 10' 


-1.9235840320 X 10' 


-1.9681933480 X 10' 


b2 


-4.9359948916 x 10' 


-2.6259765669 x 10' 


6.6356770418 x 10" 


7.1592304894 X 10" 


7.8432338862 x 10" 


bi 


-2.0507457345 x 10' 


-3.3913231244 x 10" 


1.3837014244 x 10' 


1.3979846789 x 10' 


1.4343371719 x 10' 


b4 


3.6291414100 x 10' 


2.5464049259 x 10' 


7.8618686503 x 10" 


7.3730033289 x 10" 


6.8679571684 x 10" 




4 371 3955958 x 10' 


1 6166653183 x 10' 


-8 0358407914 x 10" 


-8 653201 5972 x 10" 


-9 6695240617 x 10" 


bi 


-5.2721158916 x 10' 


-2.8335792233 x 10' 


-4.8006902023 x 10" 


-3.6779369765 x 10" 


-2.2948815920 x 10" 


h 


1.4831730999 x 10' 


8.6649789093 x 10" 


2.8424870765 x 10" 


2.4362978935 x 10" 


2.0390764869 x 10" 


b. 


-1.3230829528 x 10" 


-7.9156543417 x 10"' 


-3.2207517000 X 10"' 


-2.7645228408 x lO"' 


-2.3641043114 X 10"' 




1.7474071394 x 10"' 


1.0596166967 x 10"' 


4.6227764509 x 10^^ 


3.9732841793 x 10^^ 


3.4157538873 x 10"^ 




5.0793294231 x 10' 


5.9525977879 x 10' 


6.8805689403 x 10' 


7.8512293219 x 10' 


8.8558644431 x 10' 




12s 


13s 


14s 


15s 


16s 




-3.9188326323 x 10' 


-3.9670102478 x 10' 


-4.0057506143 x 10' 


-4.0380409686 x 10' 


-4.0661454371 x 10' 




-3.5484078761 x 10"' 


-4.0844617686 x 10"' 


-4.5428570381 x 10"' 


-4.6335483910 x 10"' 


-4.9050805652 x 10"' 


bo 


4.9431541043 x 10' 


2.9197464029 x 10' 


1.7971514946 x 10' 


1.5982980641 x 10' 


1.1724743879 x 10' 


bi 


-1.9750664049 x 10' 


-1.9723076238 x 10' 


-1.9651279272 x 10' 


-2.0606406024 x 10' 


-2.1081283167 x 10' 


b2 


8.1838767684 x 10" 


8.3673624537 x 10" 


8.4788048970 x 10" 


9.3338234414 x 10" 


9.8388495127 x 10" 


b. 


1.4326519992 x 10' 


1.4291867041 x 10' 


1.4206704409 x 10' 


1.4985604253 x 10' 


1.5395638329 x 10' 


b4 


6.4270125576 x 10" 


6.0967371906 x 10" 


5.8366404588 x 10" 


5.7371 156374 X 10" 


5.5633379574 x 10" 


he 


-1 0032296167 x 10' 


-1 0363858125 x 10' 


-1 051 8387038 x 10' 


-1 1546724399 x 10' 


-1 2242342508 x 10' 

1 • ^^^^^^^^ \_/L) /\ L yj 


bi 


-1.4087194335 x 10" 


-6.8184121010 X 10^' 


-2.9722754662 x 10"' 


1.1575092351 x 10^' 


4.8942663430 x 10"' 


h 


1.8482968593 x 10" 


1.8529967851 x 10" 


2.1519310722 x 10" 


2.2891230659 x 10" 


2.6782065008 x 10" 


b. 


-2.1185600113 X 10-' 


-1.9498101086 X 10"' 


-1.8067549402 x 10"' 


-2.0119354736 X 10"' 


-1.8223251218 x 10"' 


b9 


3.0555243512 x 10"^ 


2.7915669075 x lO"- 


2.5650331467 x 10"^ 


3.2820939449 x 10"^ 


3.0707648484 x 10"- 




9.8879969305 x 10' 


1.0942707459 x 10" 


1.2016192323 x 10^ 


1.3105436054 X 10^ 


1.4208118814 X 10" 




17s 


18s 


19s 


20s 




ao 


-4.0917355907 x 10' 


-4.1160166944 X 10' 


-4.1398378313 x 10' 


-4.1637819958 X 10' 






-4.7538949275 x 10"' 


-4.4937224788 x 10"' 


-4.4703263317 x 10"' 


-4.8244405309 x 10"' 




bo 


1.2547183851 x 10' 


1.4317441971 x 10' 


1.3339226887 x 10' 


9.1516183240 X 10" 




bi 


-2.1972421315 x 10' 


-2.2875778234 x 10' 


-2.3839004652 x 10' 


-2.3934619536 x 10' 




b2 


1.0636513367 x 10' 


1.1277442747 x 10' 


1.2219984677 x 10' 


1.2479746639 x 10' 




b. 


1.5943304131 x 10' 


1.6640858983 x 10' 


1.7332324934 x 10' 


1.7384969031 x 10' 




bA 


5.4874376556 x 10" 


5.6104613925 x 10" 


5.3841616235 x 10" 


5.1298737644 x 10" 




bs 


-1.2655779187 x 10' 


-1.3260181912 X 10' 


-1.4146776145 x 10' 


-1.4462821999 x 10' 




bi 


4.9579530934 x 10"' 


5.2627269447 x 10"' 


1.0461078949 x 10" 


1.3716208664 x 10" 




h 


2.5916598871 x 10" 


2.4342837775 x 10" 


2.3937323753 x 10" 


2.8470359883 x 10" 




bi 


-2.2938505101 x 10"' 


-2.7807376757 x 10"' 


-2.8519500457 x 10"' 


-2.4816096711 x 10^' 




h 


4.2155577449 x 10"^ 


5.4040208750 x lO"- 


5.6731652666 X 10"- 


4.9205565009 x 10"^ 




r(l/2) 


1.5322255378 x 10^ 


1.6446312975 x 10^ 


1.7579037822 X 10^ 


1.8719377215 x 10^ 





20 Chluba and Sunyaev: Two-photon transitions in hydrogen and cosmological recombination 



Table B.2. Non-linear fitting coefficients for the non-resonant nd — » Is emission spectra within the frequency range 10 < y < 
0.999 for n < 20. 





3d 


4d 


5d 


6d 


7d 




-1.0484392971 x 10' 


-1.4527047719 x 10' 


-1.7350074013 x 10' 


-1.9527421713 x 10' 


-2.1271595034 x 10' 




2.6115596268 x lO"' 


8.1050542323 x 10"' 


4.6206291704 x 10^' 


2.3419462081 x 10^' 


6.7999906938 x 10"^ 


bo 


5.3546036998 x 10' 


1.7065087104 x 10^ 


2.0902445563 x 10^ 


1.4535416639 x 10' 


1.1921878150 X 10^ 


bi 


-9.7319568527 x 10' 


-4.2936500531 x 10^ 


5.8562543507 x 10' 


5.9612688683 x 10' 


3.0991642200 x 10' 


b2 


1.0624929126 x 10^ 


7.1072286522 x 10^ 


-7.8501054392 x 10' 


-8.0969461685 x 10' 


-5.2470184998 x 10' 


bi 


6.0569007312 x 10' 


-3.6472867939 x 10' 


-5.8200368181 x 10' 


-5.3437128442 x 10' 


-2.3362093646 x 10' 


b4 


-4.1091307521 x 10' 


-4.4275350080 x 10" 


4.8846690629 x 10' 


4.9735078265 x 10' 


3.7636220223 x 10' 




-8 3346797794 x 10' 


7 4703670467 x 1 0^ 


1 1981858951 x 10^ 


9 824784491 5 x10' 


4 8776486304 x 10' 


bi 


6.9608223453 x 10' 


-6.6433743910 x 10' 


-1.1364879196 X 10^ 


-9.8596391243 x 10' 


-5.6611227459 X 10' 


h 


-1.7046798510 x 10' 


1.8959198585 x 10' 


3.1841053645 x 10' 


2.6560347162 x 10' 


1.5852011862 x 10' 


b. 


1.6052249571 x 10" 


-3.1 192847294 X 10" 


-3.0563352329 x 10" 


-2.3898992832 x 10" 


- 1.4245686402 X 10" 


b9 


-2.2422235582 x lO"' 


5.0786391006 x 10"' 


4.1302480839 x 10^' 


3.1693775497 x 10^' 


1.8991428606 x 10"' 




2.0915177571 x 10' 


3.4902680010 x 10' 


4.4901877444 x 10' 


5.2879363207 x 10' 


5.9763238933 x 10' 




8d 


9d 


lOd 


1 Id 


12d 


ao 


-2.2687581599 x 10' 


-2.3840999089 x 10' 


-2.4781138994 X 10' 


-2.5548956191 x 10' 


-2.6179506753 x 10' 




-4.0416246997 x lO'^ 


-1.4739834639 X 10"' 


-2.3472013799 x 10"' 


-3.0648469505 x 10"' 


-3.6629270128 x 10"' 


bo 


1.8805941117 X 10^ 


1.0632962733 x 10^ 


6.0922960536 x 10' 


3.6160368635 x 10' 


2.2144091289 x 10' 


bi 


-1.8483714555 x 10' 


-1.8910499138 x 10' 


-1.9083232938 x 10' 


-1.9357145527 x 10' 


-1.9508239494 x 10' 


b2 


5.8156915967 x 10" 


6.5026790075 x 10« 


6.8739045830 x 10" 


7.3456517995 x 10" 


7.6357892226 x 10" 


bi 


1.3504738304 x 10' 


1.3795221365 x 10' 


1.3924013399 x 10' 


1.4140788681 x 10' 


1.4314213871 x 10' 


b4 


8.5177417832 x 10" 


8.0023959097 x 10" 


7.6277097156 x 10" 


7.2310576094 X 10" 


6.9389854602 x 10" 




-7 0817594819 x 10" 


-7 91 16375657 x 10'' 


-8 4667321 100 x 10" 


-9 1736549163 x 10" 


-9 8020559467 x 10" 


bi 


-6.3232881999 x 10" 


-5.0886630492 x 10" 


-4.1097110628 x 10" 


-3.0436236228 x 10" 


-2.1203514774 x 10" 


h 


3.4340616701 x 10" 


2.9768149183 x 10" 


2.6171546155 x 10" 


2.2949189350 X 10° 


2.0887633671 x 10" 


b. 


-3.9664209528 x lO"' 


-3.4089403622 x 10"' 


-2.9766246686 x 10"' 


-2.6063687833 x 10"' 


-2.3211600328 x 10"' 




5.7497854559 x lO^^ 


4.9213789949 x 10^- 


4.2727779792 x 10"- 


3.7171055385 x lO^^ 


3.2750870685 x 10"^ 




6.6011443629 x 10' 


7.1866675229 x 10' 


7.7466332091 x 10' 


8.2892729834 x 10' 


8.8197557986 x 10' 




13d 


14d 


15d 


16d 


17d 


ao 


-2.6702542838 x 10' 


-2.7142749615 x 10' 


-2.7520108750 x 10' 


-2.7850443691 x 10' 


-2.8146062687 x 10' 


P 


-4.1655350067 x lO"' 


-4.5817894101 x 10"' 


-4.9230994423 x 10"' 


-4.8799319840 X 10"' 


-5.1046103153 x 10"' 


bo 


1.4031825050 x 10' 


9.3536034158 x 10" 


6.5800387184 x 10" 


6.5003185806 X 10" 


5.0269776925 x 10" 


bi 


-1.9675880116 x 10' 


-1.9851518345 x 10' 


-2.0183574589 x 10' 


-2.0542531107 x 10' 


-2.0953818322 x 10' 


b2 


8.0418685007 x 10" 


8.3266612119 x 10" 


8.8233502592 x 10" 


8.8786234533 x 10" 


9.3483218939 x 10" 


b. 


1.4344987353 x 10' 


1.4507266322 x 10' 


1.4666068986 x 10' 


1.5088704084 X 10' 


1.5370738183 x 10' 


b4 


6.5426767345 x 10" 


6.3117201027 x 10" 


6.0075843919 x 10" 


6.3138200610 X 10" 


6.1350437791 x 10" 


he 


-1 0108197620 X 10' 


-1 05871001 58 X 10' 


-1 0920565892 x 10' 


-1 1277066540 x 10' 


-1 17144Q7606 X 10' 


bi 


-1.4725216928 x 10" 


-9.1174802515 x 10^' 


-6.0110879243 x 10"' 


-8.1971758897 x 10^' 


-6.0125221865 x 10"' 


h 


2.0501069105 x 10" 


2.1660824963 x 10" 


2.4592074550 x 10" 


2.6330961336 x 10" 


2.9489602434 x 10" 


b. 


-2.0980025116 X 10"' 


-1.8682379725 x 10"' 


-1.5933095829 X lO"' 


-2.0505355745 x 10"' 


-1.7710483631 x 10"' 


b9 


2.9129276734 x 10-^ 


2.5589737997 x 10"^ 


2.1812553698 x lO"- 


3.3518563369 x lO^^ 


2.9819941280 x 10"- 




9.3414593519 X 10' 


9.8566665797 x 10' 


1.0366966228 x 10^ 


1.0873478648 x 10^ 


1.1377060842 x 10" 




18d 


19d 


20d 






ao 


-2.8416405767 x 10' 


-2.8668633507 x 10' 


-2.8908128993 x 10' 






P 


-5.2741247590 x 10"' 


-4.8873340484 x 10"' 


-4.7794020321 x 10"' 






bo 


4.1457341691 x 10" 


5.2856466775 x 10" 


5.2735742805 x 10" 






bi 


-2.1898431473 x 10' 


-2.2963670035 x 10' 


-2.3833638246 x 10' 






b2 


1.0391335216 x 10' 


1.1284085887 x 10' 


1.2121030366 x 10' 






b. 


1.6037472635 x 10' 


1.6719978965 x 10' 


1.7273259137 x 10' 






b4 


5.7762172996 x 10" 


5.7383329948 x 10" 


5.5675791948 x 10" 






bs 


-1.2705538272 x 10' 


-1.3194456779 x 10' 


-1.3752468820 x 10' 






bi 


-1.8363323426 x lO"' 


7.2380465060 x 10"^ 


3.0584919862 x 10"' 






h 


3.1624501247 x 10" 


2.7857591072 x 10" 


2.7454079863 x 10" 






bi 


-1.4050292467 x lO"' 


-2.4398422587 x 10"' 


-2.7060311659 x 10"' 








2.5244597061 x lO"- 


4.6094366553 x 10"^ 


5.2504979483 x 10"^ 






rw2) 


1.1878311140 X 10- 


1.2377675855 x 10" 


1.2875544820 x 10^ 







